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Using numerical simulations we show that, in the presence of random quenched disorder, sliding 
superconducting vortices and Wigner crystals pass through a variety of dynamical phases when an 
additional transverse driving force is applied. If the disorder is weak, the driven particles form a 
moving lattice and the transverse response shows a devil's staircase structure as the net driving 
force vector locks with the symmetry directions of the moving lattice, in agreement with the pre- 
dictions of Le Doussal and Giamarchi. For strong disorder, and particularly for smoothly varying 
potential landscapes, the transverse response consists of a sequence of disordering transitions with 
an intervening formation of stable channel structures. 

PACS numbers: 74.25.Qt 



I. INTRODUCTION 

A number of systems, including moving vortex lat- 
tices in type-II superconductors with disordered pinning 
P) 0) @1' sliding charge-density waves Q, and driven 
Wigner crystals [5], can be modeled as moving lattices 
interacting with random quenched disorder. If the exter- 
nal driving force is weak the particles are pinned, while 
at high external drives the particles can slide. There is 
extensive evidence that in the sliding phase the parti- 
cles organize into a moving anisotropic crystal or moving 
smectic phase and that the motion is confined to well de- 
fined channels [U, i, i, i, 0, B 111 • A remarkable property 
of these channels is that if an additional external drive 
is applied in the transverse direction the moving^system 
exhibits a finite transverse depinning threshold [2] . This 
transverse depinning threshold has been confirmed in a 
number of simulations [H, 0, S and experiments 

A transverse depinning barrier is also present for parti- 
cles moving over a periodic substrate. In this type of sys- 
tem, the transverse velocity force curves exhibit a devil's 
staircase structure as the particle motion locks to the 
symmetry directions of the periodic substrate pdl. fl^. [l3| . 
This effect occurs even for a sing le particle moving over 
a periodic substrate [13, [H, In contrast, a single 

particle moving over a random substrate will not show a 
transverse depinning threshold or a devil's staircase re- 
sponse. 

In Ref . [H , a finite transverse depinning force was pre- 
dicted for a vortex lattice, and it was also conjectured 
that a moving lattice would align with the direction of 
the initial longitudinal drive in order to minimize power 
dissipation. If the moving lattice remains aligned in this 
initial orientation as a transverse drive is applied, the 
transverse velocity force curves could exhibit a devil's 
staircase response as the applied force vector locks with 
the triangular symmetry of the moving lattice. In con- 
trast to a system with a periodic substrate, here the sym- 
metry breaking which leads to the preferred directions of 
motion comes from the moving lattice itself. Numerically 



observing the devil's staircase in a system with random 
disorder can be difficult since large systems and long time 
averages of the transverse velocity are needed to reduce 
the fluctuations. An added difficultly with observing this 
type of devil's staircase is that fine increments of the ap- 
plied transverse force are required to resolve the higher 
order steps. 

Another possibility for the dynamics at the transverse 
depinning transition which has not been considered is 
that the moving lattice structure could become strongly 
disordered at the depinning transition, leading to plastic 
distortions. If such a plastic transverse depinning occurs, 
it is natural to ask if the particles would reorder into a 
tilted channel structure, causing a new transverse barrier 
to develop. 



II. SIMULATION 

A. Vortex System 

To address these issues we consider two types of two- 
dimensional (2D) systems with periodic boundary condi- 
tions. The first is a vortex system containing 7V„ vortices 
and Np short range random pinning sites. The vortices 
are modeled as point particles and their dynamics evolve 
via an overdamped equation of motion, given for vortex 
i as 

7?^ = Fr+Fr+F,, (1) 

where rj is the damping constant. The vortex- vortex 
interaction force [13, [11 is F[ = Yj^=^j foKi{nj / X)T^j, 
where Ki{r) is the modified Bessel function, rij = Ir— rj| 
is the distance between vortices i and j located at 
and Tj, Yij = (r^ — rj)/rij, the unit of force is /o = 
<I'§/27r/iO'^'^) ^i'o is the flux quantum, and the London pen- 
etration depth A is the unit of length. This interaction is 
appropriate for stiff 3D vortices in relatively thick films 
where the film thickness d is much larger than A but is 
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smaller than the Larkin-Ovchinnikov correlation length 
[T9I [20j . Time is measured in units of t = 77//0 and a 
molecular dynamics (MD) time step is equal to 0.02t. We 
find identical results if smaller time steps are used. For 
a O.lfim thick crystal of NbSeo^^/o = 6.78 x lO^^N/m 
and ?7 = 2.36 x 10"" Ns/m j2l|, so that t = 0.35/iS. 
The pinning force arises from pinning sites which are 
modeled as attractive parabolic potentials of radius Vp 
and strength /pi Ff = - J2k=i Ip{fik/rp)Q{rp - rik)Y^k, 
where rik — \Yi — Vk\ is the distance between vortex i 
and pin k located at r^, Tik = (r^ — rk)/rik, and Q is 
the Heaviside step function. The pinning interaction is 
intended to represent columnar pinning [2^ . The net 
external driving force F^; is comprised of the longitudi- 
nal force applied in the x-direction, F^ = ^it^' ^^'^ 
transverse force applied in the y-direction, F^*" = Fj'^y. 
The vortex (pin) density is given by py {pp). The sys- 
tem is prepared by simulated annealing, after which a 
longitudinal force is gradually applied in increments of 
1 X 10~'*/o up to a specific value. The final longitudinal 
driving force is held fixed and then the transverse drive is 
applied in small increments. We measure the transverse 

velocity V"^^ = {N~^ J2i^=i '^i ' y)' where is the veloc- 
ity of vortex i. We also measure the fraction of sixfold 
coordinated vortices, Pg = -^iT^ '^(6 ~ Zi), where 

Zi is the coordination number of vortex i as determined 
from a Voronoi construction. 



B. Sliding Wigner Crystals 

The second system we consider is sliding Wigner crys- 
tals composed of electrons. Here we use the same 
Langevin type simulations which were employed previ- 
ously to identify the sliding states in this system . The 
electron-electron interaction force F[ = —VUl, where 

ui^y — - (2) 

and ri(j) is the location of election The electron 

crystal has lattice constant cq. The disorder comes from 
Np positively charged impurities placed at a distance d = 
0.9ao out of plane with the form Ff — — VC/f , where 

2 

£/f = -V^ (3) 

iv) 

Here is the in-plane projection of the location of pin 
J. The number of electrons Ne equals the number of im- 
purities Np. Both the electron-electron and electron im- 
purity interactions are long rang e, and we evaluate them 
with a fast-converging sum [23[. To measure the trans- 
verse response, we follow a similar procedure as in the 
vortex system. 



III. TRANSVERSE DEPINNING FOR VORTEX 
SYSTEMS 

A. Definition of Pinning Regimes 

We first consider a vortex system for a sample of size 
L = 24A with p„ = = O.TS/A^, Vp = 0.2A, and 
fixed longitudinal drive = 0.012/o. For this value 
of and pinning density we find three regimes of the 
transverse response for varied pinning strength. For 
fp < O.O85/0 the system is in what we term the weak 
pinning limit, and the equilibrium state at F^ = 0.0 is 
free of dislocations. At Fjj^ = 0.012/o, there is a finite 
transverse depinning threshold and the transverse depin- 
ning occurs without the generation of additional defects. 
In the intermediate regime 0.085 < fp/ fo < 0.13, the 
equilibrium Fj' — state contains dislocations, and for 
F^ — 0.012/0 the transverse depinning is plastic and is 
accompanied by the generation of defects. Under a lon- 
gitudinal drive, for 0.085 < fp/fo < 0.106 the vortices 
dynamically reorder into a defect-free lattice, while for 
0.106 < fp/fo < 0.13 the vortices undergo a partial dy- 
namical reordering into a moving smectic state contain- 
ing a finite number of dislocations that are all aligned 
with the driving direction, as seen in previous simula- 
tions d, 0] ■ The number of defects in the moving smectic 
state gradually increases with increasing fp until each row 
of vortices is completely decoupled from the neighboring 
rows. A weak smectic state contains only a few disloca- 
tions so that several adjacent rows are coupled but rows 
separated by longer length scales are decoupled. In the 
strong pinning regime fp > 0.13/o, the vortices flow plas- 
tically through the system and the moving channel struc- 
ture is destroyed so that the dislocations are no longer 
aligned with the drive. In this case there is no transverse 
depinning barrier. We have tested system sizes up to 
L = 92A and we find only dislocation free lattices in the 
equilibrium state for fp < 0.085/o. It is possible that for 
much larger systems than we can simulate, dislocations 
will eventually appear in the weak pinning systems. 



B. Elastic Transverse Depinning 

We consider a sample with fp = O.O6/0 in the weak pin- 
ning regime where the equilibrium lattice contains no dis- 
locations. We fix the longitudinal drive to Fj' — 0.012/o 
and increase Fj'^ from zero in increments of 4 x 10~^/o 
every 1.5 x 10^ MD time steps. In order to resolve the 
steps in the transverse response, we average the trans- 
verse velocity V'^'^ over 9 x 10'' MD steps. The appear- 
ance of the curve does not change if V'^'^ is averaged 
over 1.5 x lO'' MD steps instead. A typical simulation 
runs for a total of more than 5 x 10^ MD steps. When 
the net driving force vector aligns with a symmetry di- 
rection of the vortex lattice, we find that the transverse 
response locks to that direction over a range of Fj^ . For 
a moving triangular lattice, this locking occurs whenever 
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FIG. 1: (a) The transverse velocity V'^^ vs the transverse 
applied force Fj^ for a vortex lattice moving over random 
disorder of strength fp = O.O6/0 with fixed longitudinal driv- 



ing force 



0.012/0. (b) dV'^''/dFj'' vs Fj'' for the same 



system, (c) V"^"^ vs FJ^ showing the evolution of the (1,2) 
step for /p = 0.02, 0.04, 0.06, and O.O8/0, from top to bottom, 
(d) Step widths Ws in Fj'' vs /p for the (0,0), (1,3), (1,2), 
(3,4), (3,3) and (4,4) steps, from bottom to top. 



FJ"^ /Fj' ~ \/3m/ (2n+ 1) where n and m are integers. In 
such a devil's staircase, the most pronounced steps occur 
when both n and m have low values. 

The plot of V'^'' vs Fj'' in Fig. [^a) shows that a 
transverse critical depinning threshold F^^' occurs at 
FJ^ = O.OOO65/0. It is important to note that the vor- 
tex lattice remains aligned in the a;-direction during the 
entire simulation. For increasing Fj"^ , a number of steps 
appear in V'^^ corresponding to the locking of the net 
applied force vector with a symmetry direction of the 
lattice. A prominent step occurs near Fj^ = 0.0043/o 
which corresponds to the case of (m, n) = (1, 2). The step 
at FJ"" — 0.003/0 corresponds to the (1,3) state, and the 
(2,3) state appears at Fj^ = 0.007/o. There are also a 
number of smaller steps at higher values of n and m, con- 
firming the existence of a devil's staircase response. In 
Fig. [Tfb), dV'^^ / dFj^ vs FJ^ more succinctly highlights 
the steps. The peaks in Fig. 1(b) correspond to points at 



which the curvature in V^"^ vs Fj"^ increases just before 
and after a plateau region. The most pronounced steps 
occur for (m, n) where the value of both m and n are 
small, which is typical for a devil's staircase structure. 
This result confirms the prediction of Giamarchi and Le 
Doussal that if the moving lattice stays aligned with 
the original applied longitudinal drive, a devil's staircase 
response occurs due to the matching of the symmetry 
directions of the lattice with the applied drive. 

The width of a given step Ws depends on the pinning 
strength and density. In Fig. [Tfc) we illustrate the evo- 
lution of the (1, 2) step for fp/ fo = 0.02, 0.04, 0.06, and 
0.08, showing the growth of Wg with fp. We also find 
that the step widths increase when the pinning density 



FIG. 2: (a) V ^ vs F^^ for a moving vortex lattice with the 
same parameters as in Fig. 1(a), but with fp = O.IO6/0. (b) 
The fraction of six-fold coordinated vortices Pe vs FJ'' for the 
same sample showing the correlations of the drops in Pg with 
the transitions in V'^'' . 



Pp is increased. In Fig. [TJd) we plot the evolution of 
Ws as a function of pinning strength for {m,n) = (0,0), 
(1,3), (1, 2) (3, 4), (3, 3) and (4, 4). We have verified that 
these results are robust up to systems of size L = 60A 
provided that the lattice remains free of dislocations. In 
two-dimensional pinned systems it has been argued that 
dislocations eventually appear at large distances even for 
weak disorder. These dislocations would limit the reso- 
lution of the devil's staircase, which should remain ob- 
servable as long as the distance between dislocations is 
comparable to or greater than the system size. 



C. Plastic Transverse Depinning 

In the upper range of the intermediate pinning strength 
regime 0.106 < fp/ fo < 0.13, the longitudinally driven 
vortices form a moving smectic state and the trans- 
verse depinning transition is plastic. Here the transverse 
velocity response is characterized by a series of sharp 
disordering-reordering transitions, giving rise to consec- 
utive pronounced steps in V'^^ , as shown in Fig.[2la) for 
a system with fp — O.IO6/0. This ordering-disordering 
velocity-force response, which is distinct from a devil's 
staircase structure, occurs when the pinning is strong 
enough that the particles move in channels that are 
separated by a number of aligned dislocations, which 
is indicative of a smectic ordering. We also find two 
prominent peaks in the structure factor ^(k) consis- 
tent with a smectic state, as illustrated in Fig. [3] Here 
S{k) = iV-i|X;^"iexp(-ik • r,)|2. In addition to the 
two prominent peaks. Fig. [3] shows the presence of four 
smaller peaks. Since fp = O.IO6/0 falls at the edge of the 
regime where the moving smectic state forms, there are 



FIG. 3: Structure factor S'(k) for a system with fp = O.IO6/0 
in the moving smectic state at FJ^ = 4 x 10~^/o and — 
0.012/0. The central peak is not shown. 



only a small number of dislocations present in the mov- 
ing smectic state and correlations between neighboring 
channels of vortices occur over a width of several chan- 
nels. These correlations produce the four smaller peaks 
in S(k). For larger system sizes L these four small peaks 
shrink in magnitude while the two prominent peaks re- 
main prominent. In Ref. the prominent peaks were 
observed to decay algebraically with system size and the 
state was termed a moving transverse Bragg glass. It is 
beyond the scope of this work to confirm the existence 
of a moving transverse Bragg glass. For larger values 
of fp, the dislocation density in the moving smectic state 
increases, the correlation length between the moving vor- 
tex channels decreases to the width of a single channel, 
and only two peaks appear in S'(k). For fp > 0.13/o 
the channel structure is destroyed and the vortices move 
plastically. 

Since the transverse depinning process at fp — O.IO6/0 
is plastic, a large number of transient dislocations are 
induced at the depinning transition which destroy the 
well ordered channel structures. In Fig. [H^b) we plot the 
fraction of sixfold-coordinated vortices Pg for the same 
system in Fig. [Ifa), indicating that the changes in the 
behavior of F^'' coincide with changes in Pg- In the 
regions where V'^^ is increasing linearly with drive, Pg 
drops and shows rapid fluctuations. Here, the ID chan- 
nel structure of the moving lattice is destroyed. We note 
that for 0.085 < fp/fo < 0.106 we do not observe any 
dislocations in the longitudinal moving lattice; however, 
since the system is close to the value of fp at which dis- 
locations appear in the longitudinally moving lattice, at 
the transverse depinning transition the extra strain on 
the lattice is enough to induce some dislocations and the 
system also depins plastically. 




(e) 

FIG. 4: (a,c,e): Voronoi construction for the vortex config- 
urations for the system in Fig. [2] taken from a single snap- 
shot. Small dots: sixfold-coordinated vortices; large dots: 
nonsixfold-coordinated vortices. (b,d,f): Vortex trajectories 
(black lines) during a fixed time span. (a,b): Transversely 
pinned state at FJ"^ = 4 x 10~*/o where a clear channel struc- 
ture appears. (c,d): The transition state at FJ^ — O.OOll/o. 
(e,f): The state at the first plateau in V'^'' at FJ'' = 0.00125/o 
where the system has reorganized into a tilted channel struc- 
ture. 



To further illustrate the change in lattice structure that 
occurs on the steps in the intermediate pinning regime, in 
Fig. [31^ a) we show the Voronoi construction of the vortex 
lattice from the system in Fig. [2] for a transversely pinned 
state at Pj'' = 4 x IO^^/q. There are a small number of 
dislocations present, indicated by large dark circles, and 
all of the dislocation pairs are aligned with the direction 
of the drive (to the right in the figure). The vortex lattice 
is also aligned with the driving direction. In Fig. |4{b) 
the vortex positions (black dots) and trajectories (black 
lines) indicate that the motion is confined to well defined 
ID channels. Above the transverse depinning transition, 
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illustrated in Fig. [4ljc,d) for Fj"^ = O.OOll/o, numerous 
dislocations are present as indicated by Fig. ID^c), while 
the trajectories in Fig. \M.d) show that the channeling 
effect has been lost. This disordered type of flow appears 
both at the transitions among the higher steps as well as 
in the regions where V'^^ is increasing linearly. 

In Fig. Elja), the velocity jumps to a finite value and 
becomes locked onto a plateau for O.OOll/o < Fj^ < 
0.0017/0. Here and on the other plateaus in V'^^, the lat- 
tice moves in a channel structure which is aligned at an 
angle to the x direction, such as that shown in Fig.|He,f) 
for Fd = 0.00125/0. As FJ^ increases, the system un- 
dergoes repeated plastic depinning transitions, each of 
which appears as a drop in Pg. In the plastic flow region 
that follows each of these depinning transitions, the vor- 
tices move in the direction of the net force vector. As 
the transverse drive increases to a new plateau, the lat- 
tice reorders into a new channel structure which sets up 
a new transverse barrier to further increases in V"'"^ . 

The channels require some time to reform, so the width 
of the reordering transition as a function of transverse 
driving force is dependent on the sweep rate of Fj^ . For 
higher sweep rates, the width of the disordered regions 
between plateaus grows. In the adiabatic regime, the 
steps would be much sharper, suggesting that the transi- 
tions are first order in nature. We also find hysteresis in 
the steps when we reverse Fj^ . For very strong pinning, 
fp > 0.13/0, and fixed Fj' — 0.012/o, the channel struc- 
tures are destroyed and the vortices flow plastically for 
all transverse drives; however, for higher values of Fj' , a 
partially ordered phase appears and the transverse depin- 
ning occurs in the same manner as in Fig. [S] In general, 
for stronger disorder the longitudinal moving state con- 
tains numerous dislocations and the transverse depinning 
occurs in a step like fashion. 

The transition between devil's staircase behavior at 
weak pinning and the order-disorder transitions associ- 
ated with intermediate pinning strength is determined 
by the appearance of dislocations in the longitudinally 
driven lattice. In general, a devil's staircase appears 
when the disorder is weak enough that there are no dis- 
locations present at any drive. For this reason, no evi- 
dence of a devil's staircase appeared in the simulations 
of Ref. which showed the existence of a transverse 
depinning barrier. These simulations were performed in 
the intermediate pinning limit, as indicated by the linear 
scaling of fc with the pinning strength in Ref. [6| and 
the presence of aligned dislocations in the lattice. The 
ordering-disordering steps that we flnd in the intermedi- 
ate pinning limit were not observed in Ref. Q since the 
resolution of the transverse force versus velocity curves 
was too low. The existence of a weak pinning regime may 
depend on system size. It has been argued that for any 
2D system with random disorder, dislocations will always 
appear at sufficiently long length scales Rd even if the 
disorder is weak. In our simulations we have considered 
numerous different system sizes, and for fp < 0.085/o we 
only observe dislocation free lattices. Due to computa- 
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FIG. 5: Fj"" versus fp for = 0.012/o showing the three 
different response regimes. For weak pinning, fp < 0.085/o, 
the transverse depinning is elastic. For intermediate pinning, 
0.085 < fp/fo < 0.13, the transverse depinning is plastic. 
There is no transverse depinning barrier for strong pinning 
fp > 0.13/0. 



tional constraints we could only simulate systems of size 
L = 92A and less. It is possible that for much larger 
systems, dislocations will eventually appear. If enough 
dislocations are present to destroy the orientational or- 
der, the devil's staircase behavior would be lost when the 
well defined lattice directions of the moving lattice are 
lost. It is then possible that for system sizes L < Rd the 
devil's staircase behavior will occur, while for L > Rd, 
the strong pinning behavior will be observed. On the 
other hand, if the system is in 3D with point pinning but 
with stiff vortices, than a 3D moving Bragg glass could 
form without dislocations, permitting a devil's staircase 
behavior to be observed for all lengths. We also note 
that the argument that dislocations will appear at some 
length scale Rd was made in the context of equilibrium 
systems. In the moving vortex case it may still be pos- 
sible to obtain a defect free lattice in two dimensions as 
suggested by Koshelev and Vinokur [l|. 

In Fig. 5 we plot the transverse depinning barrier F^^' 
vs fp and indicate the three different pinning regimes 
which can be distinguished by the nature of the trans- 
verse response. In the weak pinning regime for fp/fo < 
0.085, the transverse depinning is elastic and a devil's 
staircase forms when Fj^ is swept. Within this regime, 
F^^' increases monotonically with fp. For the interme- 
diate pinning regime 0.085 < fp/fo < 0.13, the trans- 
verse depinning is plastic and F^^ is much higher than it 
was in the elastic depinning regime. For strong pinning 
fp/fo > 0.13, the vortex lattice does not reorder when 
the longitudinal drive is applied, and the transverse de- 
pinning barrier disappears. 
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FIG. 6: {a.)V'^^ vs i^J'' for a sliding Wigner crystal with 
fixed Fj- = 0.03. (b) Pe vs FJ". (c) vs FJ". (d) The 
particle trajectories for the transversely pinned sliding state at 
FJ"" = 3 X 10"*. (e) The hysteresis in the F^'' vs FJ'' curve. 
Lower curve: increasing Fj^ . Upper curve: decreasing Fj^ . 



fluctuating force so the concept of an effective tempera- 
ture does not apply. This means that a partially ordered 
moving state can persist to drives much closer to the lon- 
gitudinal depinning threshold than in the case of short 
range pinning. 

In Fig. Etd), the particle trajectories for the case of 
long range charge disorder show that the channels are 
much more meandering than in Fig.^b). In Fig.[nt^a,b,c) 
we plot V'^^, Pq, and the longitudinal velocity for 
the sliding Wigner crystal system with fixed = 0.03. 
In this case the system initially forms a moving smec- 
tic system with Pq = 0.9. A step structure appears in 
V'^^ . On the step, the particles are flowing in channels 
and Pq shows only small fluctuations. In regions where 
V'^^ is linearly increasing, Pg rapidly fluctuates around 
a lower value, indicating that the system is more dis- 
ordered. The order-disorder transitions also affect V^, 
which decreases and shows strong fluctuations in the dis- 
ordered regimes. In the ordered regimes, there are few 
fluctuations in and there is a trend for to increase 
with FJ^ . This effect is similar to loading-unloading cy- 
cles of the type found in friction systems. In Fig[BI^e) we 
illustrate the presence of hysteresis in V'^^ vs Fj^ , which 
gives evidence that the order-disorder transitions have a 
flrst order characteristic. 



SUMMARY 



IV. TRANSVERSE DEPINNING FOR SLIDING 
WIGNER CRYSTALS ON SMOOTH POTENTIAL 
LANDSCAPES 

We now consider the effect of a smoothly varying po- 
tential substrate of the type that arises in Wigner crys- 
tal systems for charge disorder. In general, with this 
type of disorder we observe much larger relative criti- 
cal transverse depinning thresholds F^"^ / F^ than with 
short range pinning. For long range disorder F^'^/F^ 
can be as high as 0.4, while for short range disorder it 
is less than 0.1 to 0.01. The effective shaking tempera- 
ture T*'', which arises from the rapidly fluctuating force 
caused by the pinning, for particles moving over short 
range random disorder at velocity V is T^'* c>c 1/V [l[. 
At high drives, the effective temperature drops with in- 
creasing and the system can reorder. At lower drives, 
V is smaller and the system melts into the plastic flow 
phase. In the smectic phase, the transverse depinning 
barrier increases as the longitudinal drive is decreased 
Q so that the maximum transverse barrier occurs just 
before the longitudinal drive is so low that the shaking 
temperature melts the lattice and Fj*" is lost. For long 
range, smoothly disordered substrates, there is no rapidly 



In summary, using numerical simulations we have con- 
flrmed the prediction of Giamarchi and Le Doussal that a 
devil's staircase transverse response occurs for an elastic 
vortex lattice moving over random pinning when the net 
force vector aligns with the symmetry directions of the 
moving lattice. The devil's staircase response appears 
in the weak pinning regime when both the equilibrium 
state and the moving state are free of dislocations and the 
transverse depinning occurs elastically without the gen- 
eration of defects. For intermediate pinning strengths, 
where the equilibrium state is dislocated and the vortices 
form a moving smectic state under a longitudinal drive, 
the transverse depinning occurs in a series of hysteretic 
order-disorder transitions where the system alternately 
becomes disordered and then reorganizes into channel 
structures that exhibit a transverse depinning threshold. 
In this case a series of pronounced steps in the trans- 
verse response appear. In the case where the quenched 
disorder is long-range and smoothly varying, the order- 
disorder transitions dominate the response to an increas- 
ing applied transverse force. These results are applicable 
to vortices and Wigner crystals moving over random dis- 
order. 

This work was supported by the U.S. Department of 
Energy under Contract No. W-7405-ENG-36. 
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